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SUMMARY 


This  report  prescribes  a  method  for  calculating  pressure 
coefficients  and  local  Mach  numbers  for  lifting  and  non- 
lifting  pointed  bodies  of  revolution  and  for  several  special 
cases  of  blunt-nosed  bodies.  The  method,  which  utilizes  hy¬ 
brid  tandem  solutions  Involving  Generalized  Newtonian  and 
Shock-Kxpansion  theories,  provides  accurate  results  for  a 
variety  of  nose  shapes  and  fineness  ratios  over  a  range  of 
supersonic/hypersonic  Mach  numbers.  The  numerical  simplicity 
of  the  method,  which  makes  it  readily  applicable  for  quick 
hand-calculational  procedures,  was  the  prime  factor  in  its 
selection  and  publication;  the  few  existing  methods  which 
yield  accurate  results  over  a  comparable  range  of  application, 
such  as  the  method  of  characteristics,  are  not  used  extensively 
because  of  the  lengthy  numerical  calculations  involved. 

The  present  method  has  been  coopered  with  exact  solutions, 
various  pertinent  theories,  and  experimental  data  where  avail¬ 
able  and  the  overall  agreement  of  the  results  is  quite  favorable. 
The  investigation  for  lifting  bodies  indicates  the  present  method 
is  applicable  for  bodies  of  revolution  at  angles  of  attack  up 
to  about  ten  degrees. 

This  report  presents  numerical  examples  showing  stepwise 
calculational  procedures  for  obtaining  pressure  coefficient  and 
local  Mach  number  distributions  along  the  meridians  of  a  body 
of  revolution  at  angle  of  attack.  In  order  to  make  the  report 
immediately  useful  to  the  engineer  desiring  such  information, 
all  of  the  necessary  tables  and  look-up  parameters  are  in¬ 
cluded  in  the  appendices. 
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I .  INTRODUCTION 

There  are  a  variety  of  Methods  available  for  the  calculation  of  pressure  dis¬ 
tributions  over  lifting  and  non-lifting  bodies  of  revolution  in  the  supersonic- 
hypersonic  speed  regime.  The  usefulness  and  range  of  applicability  of  these 
methods  have  been  widely  discussed  in  available  literature  and  are  briefly  sum¬ 
marized  herein. 

In  general  it  can  be  said  that  the  methods  having  the  most  extensive  range 
of  applicability  involve  lengthy  calculational  procedures.  Impact  theory  pro¬ 
vides  a  qulc*  and  easy-to-apply  method  for  calculating  surface  pressures  in  the 
supersonic-hypersonic  regime  il  the  surface  is  normal  or  sufficiently  oblique  to 
the  flow.  The  generalized  shock-expansion  method  for  steady  two  dimensional  flows 
is  also  quick  and  easy  to  apply.  Its  useful  range  of  application  is  for  Mach 
numbers  greater  than  2.0  and  for  values  of  the  hypersonic  similarity  parameter, 
MQd/ l ,  greater  than  about  1.2. 

The  method  presented  herein  is  a  hybrid  combination  of  modified  impact  and 
shock-expansion  theories  which  permits  the  facile  and  rapid  hand  calculation  of 
the  pressure  and  local  Mach  number  distributions  along  any  meridian  on  the  sur¬ 
face  of  a  lifting  or  non-lifting  pointed  body  of  revolution  and  many  of  the  blunt 
nosed  shapes.  This  method,  which  has  a  widu  range  of  application,  is  Intended 
for  use  in  handbook  fashion  by  engineers  needing  such  information  for  studies 
concerned  with  aerodynamic  heating,  structural  loads,  drag,  etc. 

find  tunnel  test  data  in  the  Mach  number  range  2.0  to  4.0  are  presented 
and  are  compared  with  the  theory.  These  data  show  the  pressure  distributions 
on  a  secant  ogive  nose,  a  triple-cone  ogive  combination,  a  von  Karman  nose  shape, 
and  on  hemispheres.  These  nose  sections  are  attached  to  cylindrical  afterbodies. 


SYMBOLS  AND  NOMENCLATURE 


a 


c 

C_ 


d 

GN 


K 

l 

U 

j 

M 

N 


P 

Pt 

q 

R 

R, 


r 


r 

s 

' 

V* 

SEM 

s 

V 


*,y 

x.r.t 

y 


sound  speed 

limiting  speed  due  to  expansion  into  a  vacuum 

pressure  coefficient 

specific  heat  at  constant  pressure 

specific  heat  at  constant  volume 

diameter 

Generalized  Newtonian 

hypersonic  similarity  parameter,  MQd/t 

nose  length 

Mach  number 

critical  Mach  number 

nose  caliber 

static  pressure 

total  pressure 

dynamic  pressure 

generating  radius  for  ogival  shapes 

universal  gas  constant 

nose  radius 

radius  of  focal  point 

entropy 

parameter,  entropy/universal  gas  constant 
shock-expansion  method 
segment  length 
velocity 

body-fixed  rectangular  coordinate  system 

body-fixed  cylindrical  coordinate  system 

vertical  distance  from  the  tangent  ogive  generating  point 
to  the  centerline  of  the  ogive 


J 
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_  3 

a  angle  of  attack;  angle  between  free  stream  velocity  and 

body  centerline 

Y  ratio  of  specific  heats 

A  increment 

6  angle  between  the  free  stream  velocity  and  the  tangent  to  the 

body  surface 

0  angle  between  the  body  surface  and  the  longitudinal  axis  of 

the  body 

v  Prandtl-Meyer  flow  deflection  angle 

n  radian  measure 

c  density 

0  cos~*  (1  -  2x/t )  by  definition 

♦  radial  meridian  along  body  surface 

«  Infinity 

SUBSCRIPTS 

o  free  stream  conditions 

i  flow  quantities  for  the  zero  angle  of  attack  condition 

a  flow  quantities  related  to  the  effect  of  angle  of  attack 

b  base 

c  cone 

cn  conditions  on  the  equivalent  tangent  cone  for  a  specific 

segment 

m  matching  point 

max  maximum 

N  conditions  at  nose  vertex 

n  specific  segment  of  equivalent  tangent  body 

s.o.  secant  ogive 

stag.  stagnation 

t  total 

'  primed  value;  conditions  at  the  most  forward  point  of  a 

segment  of  the  equivalent  tangent  body  (s  -  0) 


Nomenclature  used  in  Appendix  2  corresponding  to  References  26,  27,  and 
28. 

u/c,  aVc8,  Ps/Pw.  Pw/P: 

T)/p” ,  5/0" 

p0/p,  P/P,  0o/0,  Ps/p- 

Ai  ,  Bj  , 

A3  1  A3  ,  B$  ,  Bi ,  S3  ,  S3 

s  Subscripts 
1  free  stream  conditions 

s  body  surface 

w  conditions  immediately  behind  shock  wave 

—  barred  quantities;  conditions  for  zero  angle 

of  attack 

The  remaining  nomenclature  used  in  Appendix  2  are  compatible  with  values 
previously  defined. 


parameters  given  in  Reference  26 

parameters  given  in  Reference  27 

parameters  given  in  Reference  28 

coefficients  related  to  first  order 
effects  of  angle  of  attack 

coefficients  related  to  second  order 
effects  of  angle  of  attack 
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III.  ANALYSIS  AND  RESULTS 

A.  FLOW  REGIME 

The  flow  regime  considered  herein  covers  the  Mach  number  range  from 
1.5  to  approximately  8.  In  general,  the  Mach  number  range  most  seriously  con¬ 
sidered  was  from  MQ  “  2  to  6;  thus  only  a  small  portion  of  the  hypersonic  flow 
regime  (Mq>  5)  has  been  investigated.  The  hypersonic  regime,  however,  is  basi¬ 
cally  an  extension  of  the  supersonic  regime  and  if  the  fluid  retains  the  prop¬ 
erties  of  an  ideal  gas  and  the  ratio  of  specific  heats  (cp/cy)  remains  constant, 
the  theories  developed  for  the  supersonic  continuum  should  be  valid  in  the 
hypersonic  regime.  For  extremely  high  Mach  numbers,  of  course,  the  properties 
of  the  gas  may  change  radically  and  consequently  the  theory  which  was  valid  at 
the  lower  Mach  numbers  must  be  modified  accordingly. 

B.  METHODS  OF  CALCULATING  PRESSURE  DISTRIBUTIONS 
1.  Generalized  Newtonian  Theory 

The  Newtonian  impact  theory  provides  the  engineer  with  a  simple, 
rapid,  and  easy-to-apply  method  for  determining  the  local  pressure  and  Mach 
number  distributions  over  a  wide  variety  of  nose  shapes.  The  flow  concept  as 
formulated  by  Newton  assumes  that  the  shock  wave  lies  on  the  body  surface  and 
that  the  component  of  momentum  normal  to  the  body  surface  is  transferred  to  the 
body  while  the  tangential  component  remains  unchanged;  a  condition  which  is 
reached  in  the  limit  as  Mq  approaches  Infinity  and  y  approaches  1.  This  concept 
provides  a  force  or  force  coefficient  (Cp)  which  is  dependent  only  on  the  local 
slope  of  the  surface,  i.e., 

C  -  2  sin'Q 
P 

The  above  expression  for  Cp  reduces  to  that  of  the  hypersonic  small  disturbance 
theory, 

cp  ■  (r  +  1)  8s  -  2  8* 

where  the  assumption  of  very  slender  bodies  at  high  Mach  numbers  yields: 

r—+  1  as  Mq— *>  - 

and  sin  0  ■  0  for  small  angles. 

The  Newtonian  concept  of  flow  as  stated  above  neglects  the  effects 
of  centrifugal  forces  which  arise  due  to  body  curvature.  When  the  body  curvature 


i 


i 


jfr 


wrt 


is  small  in  the  free  stream  direction,  the  centrifugal  forces  in  the  thin  layef 
of  air  between  the  shock  wave  and  body  surface  should  not  affect  the  impact 
pressures  appreciably.  When  the  body  curvature  is  large,  however,  the  impact 
pressures  might  be  significantly  altered  by  the  centrifugal  forces  which  tend 
to  relieve  the  impact  pressures.  Various  authors  have  attempted  to  modify  the 
simple  impact  theory  expression  to  account  for  centrifugal  forces  and  have  don* 

M 

so  with  little  or  moderate  success.  Busemann1  presents  the  problem  in  some  ae' 
tail  with  encouraging  results,  but  the  equation  for  determining  the  local  pres' 
sure  coefficient  becomes  somewhat  more  complicated  than  that  given  by  the 
original  simple  impact  theory. 

Since  the  Newtonian  theory  offers  a  convenient,  easy-to-apply  metl»oc* 
of  determining  pressure  distributions,  further  attempts  were  made  to  improve  up01* 
the  original  Cp  equation;  these  efforts  proved  to  .^e  acceptable  only  when  used 
with  values  of  the  hypersonic  similarity  parameter  (K)  greater  than  two.  To 
this  extent,  Lees*  suggested  the  modified  form 


C  -  C 
P  P, 


sin2© 


max 


for  use  with  blunt  bodies  where  C_ 


is  derived  from  the  normal  shock  relation5* 


*max 

By  generalizing  Lees'  modified  form,  Love3 indicated  the  above  equation  could  be 
made  to  include  pointed-nosed  bodies  when  used  in  the  form: 

C  -  C  sln*e 
p  PN  sin*©.. 

B 


where  N  refers  to  conditions  at  the  nose  vertex.  For  ©N  -  90®,  of  course,  Love'® 
generalization  reverts  to  Lees'  modified  form  but  for  ©N  values  less  than  90°  , 
an  expression  is  obtained  which  exhibits  advantages  over  the  use  of  simple  impact 
theory  (Cp  -  2  sin*©).  To  illustrate  graphically,  Love's  Generalized  Newtonian 
theory  has  been  compared  with  the  Newtonian  theory  and  exact  (method  of  charact5**^ 
istics)  solutions  for  several  tangent  ogives  in  Figure  1.  The  exact  solutions 
were  obtained  from  Rossow4  and  include  the  effects  due  to  rotation.  As  shown  if 
Figure  1,  the  Generalized  Newtonian  theory  offers  excellent  agreement  with  the 
exact  solutions  up  to  approximately  an  x/l  value  equal  to  0.6;  of  equal  interest 
is  the  theory's  apparent  independence  of  K,  the  hypersonic  similarity  parameter- 
In  contrast,  the  Newtonian  theory  is  markedly  affected  by  K  as  indicated  by  the 
shaded  portion  of  the  figure  wherein  K  varies  from  0.5  to  4;  from  this,  the 
Newtonian  theory  would  appear  to  approach  the  exact  solutions  only  for  K  >  >>  1 . 


U 

Reference  numbers  will  always  appear  as  a  superscript  numeral  to  the  right  of 
the  subject  and  the  footnote  symbol  (#)  will  always  appear  as  a  super¬ 
script  to  the  left  of  the  subject. 


Fig.  1  COMPARISON  OF  NEWTONIAN  AND  GENERALIZED- NEWTONIAN  THEORY 
WITH  EXACT  SOLUTIONS  FOR  OGIVES 

These  results  indicate  the  Generalized  Newtonian  theory  would  be  highly  suitable 
for  use  in  calculating  pressure  distributions  on  ogival  shapes  if  one  further 
improvement  could  be  made,  i.e.,  refinement  of  the  theory  to  better  predict  the 
pressure  (especially  for  the  lower  K  values)  over  the  after  40  percent  of  the 
nose  length  (an  Important  region  for  drag  consideration).  In  the  interests 
of  rapid  calculational  use,  however,  the  relative  simplicity  of  the  theory  should 
be  maintained,  and  on  this  basis,  further  improvement  of  this  type  appears  re¬ 
mote  . 

2.  Shock-Expansion  Method 

The  shock-expansion  method®  is  an  extremely  simple  method  for 
determining  the  flow  fields  about  bodies  in  supersonic  flow,  but  the  range  of 
Mach  numbers  and  fineness  ratios  which  can  be  treated  by  this  method  are  limited. 
To  apply  the  method,  the  flow  conditions  at  the  nose  vertex  c  the  body  are  ob¬ 
tained  from  conical  flow  solutions  (charts  or  tables)  and  for  flow  downstream  of 
the  nose  vertex,  the  Prandtl-Meyer  flow  equations  or  tables  are  used. 
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Comparing  solutions  obtained  by  the  shock-expansion  method  with 
exact  solutions  indicates  that  accurate  pressure  distributions  over  ogival  nose 
shapes  can  be  predicted  by  the  shock-expansion  method  provided  the  hypersonic 
similarity  parameter,  K,  is  greater  than  approximate  ly  1.2  (comparisons  of  this 
method  with  other  methods  will  be  discussed  further  in  Section  D) . 

3.  Method  of  Characteristics 


The  method  of  characteristics'’  provides  accurate  solutions  to 
the  complete  differential  equations  of  flow  by  the  use  of  numerical  and  graphica 
integration.  The  accuracy  of  the  method  is  dependent  upon  the  fineness  of  the 
"mesh"  employed  and  since  the  solution  by  hand  calculation  requires  much  tedious 
work,  the  lime  required  per  solution  becomes  prohibitive  for  most  generalized 
design  studies. 

The  method  of  characteristics  has  been  applied  to  many  basic 
types  of  flow  fields  including  many  of  those  associated  with  the  configurations 
contained  herein.  As  a  result,  this  method  has  been  used  as  a  basis  for  com¬ 
parison  when  either  experimental  data  or  pertinent  theoretical  solutions  are 
unavailable.  The  validity  of  using  pressure  distributions  from  characteristic 
solutions  as  standards  has  been  well  established  by  its  correlation  with  avail¬ 
able  experimental  data. 
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4.  Tangent-Cone  Method  — 

The  tangent-cone  method®  estimates  pressures  on  curved  bodies  by 
using  the  flow  solutions  for  cones  which  have  the  same  slopes  as  those  of  the 
body  surface  at  the  points  in  question.  Two  applications  of  this  method  are 
available  wherein  the  relative  accuracy  of  either  application  depends  on  the 
Mach  number  and  fineness  ratio  of  the  nose.  One  application  of  the  tangent-cone 
method  is  to  merely  use  the  pressure  coefficient  for  cones  whose  semi-vertex 
angles  correspond  to  those  of  the  body  surface.  This  assumes  a  different  total- 
head  pressure  for  each  segment  and,  in  addition,  cannot  predict  negative  pres¬ 
sures.  This  application  is  referred  to  as  the  tangent-cone  method  with  local 
total-head  ratio. 


The  second  application  of  the  tangent-cone  method 
vertex  total-head  ratio  across  the  bow  shock  wave  of  the  body  in 
of  the  surface  pressures  where  the  local  Mach  numbers  correspond 
cones  tangent  to  the  body  surface.  This  application  is  referred 
tangent-cone  method  with  vertex  total-head  ratio. 


uses  the  nose 
the  calculation 
to  those  for 
to  as  the 


The  tangent-cone  method  with  local  total-head  ratio  has  been  used 
briefly  herein  for  comparison  purposes  and  this  particular  application  was 
selected  since  it  offers  better  agreement  with  exact  solutions. 
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5.  Other  Methods 

There  are  various  other  methods  or  theories  available  for  calcula¬ 
ting  pressure  distributions  such  as  Linearized  theory’,  First  Order  theory8, 
Second  Order  theory8,  Slender  Body  theory1 15  ’ 1 1  ,  etc.,  each  of  which,  in  general, 
is  applicable  for  a  limited  range  of  Mach  numbers  and  nose  fineness  ratios. 
These  methods  will  not  be  discussed  since  for  the  most  part,  they  represent 
solutions  which  are  not  readily  adapted  to  quick  hand  calculational  procedures. 

6,  Present  Method 

Various  authors  in  dealing  with  high  supersonic  or  hypersonic 
flow  over  hemispherical  noses  have  successfully  used  a  combination  impact-shock 
expansion  theory  to  predict  pressure  distributions.  This  blunt  body  method  is 
discussed  in  Section  E.  The  present  method  extends  this  hybrid  type  of  approach 
to  "pointed  bodies"  in  that  Generalized  Newtonian  theory  is  used  to  obtain  the 
pressure  distribution  on  the  forward  portion  of  the  nose  and  the  shock-expansion 
method  is  used  to  obtain  the  pressure  distribution  on  the  after  section  of  the 
nose.  The  pressure  distribution  on  the  cylindrical  afterbody  is  obtained  by 
means  of  an  exponential  decay  law  derived  by  Fenter1 8  from  a  second  order  shock- 
expansion  approximation.  In  the  present  analysis,  a  method  has  been  developed 
for  the  tandem  blending  of  these  solutions  to  provide  a  rapid  and  accurate  hand 
calculational  procedure  for  obtaining  the  pressure  and  local  Mach  number  dis¬ 
tributions  over  a  variety  of  nose  curvatures.  This  method  will  be  developed  in 
a  stepwise  procedure  in  Part  1  of  Section  D  which  describes  its  application  to 
tangent  ogives. 

C.  METHOD  OF  CORRELATING  PRESSURE  DISTRIBUTIONS 


The  pressure  distributions  calculated  by  the  "Present  Method"  for 
pointed  bodies  of  revolution  are  correlated  on  the  basis  of  the  Hypersonic 
Similarity  Parameter,  K,  which  represents  the  ratio  of  free  stream  Mach  number 
to  nose  fineness  ratio.  The  hypersonic  similarity  rule  states  that  for  related, 
pointed,  axially  symmetric  bodies  of  revolution,  the  pressure  distributions  (in 
terms  of  (p  -  PQ)/P0  or  the  normalizing  parameter  C  /C  )  depend  only  on  the 

similarity  parameter  K.  Thus,  if  the  pressure  distribution  for  a  given  body  is 
known,  the  pressure  distributions  for  geometrically  similar  bodies  are  identical 
provided  the  two  bodies  have  the  same  value  of  the  similarity  parameter,  K. 
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The  range  of  application  of  this  rule  for  ogives4  is  shown  in  Figure  2  and  as  noted, 
the  rule  is  normally  not  valid  for  Mach  numbers  or  fineness  ratios  less  than 
approximately  2.  The  present  method  illustrates  this  limit  may  be  relaxed  slightly 
to  include  a  nose  fineness  ratio  of  1.5  for  K  ■  2. 

The  hypersonic  similarity  ~*ule  has  been  successfully  applied  to  sever¬ 
al  special  cases  included  herein  which  theoretically  have  blunt  noses  (infinite 
slope  at  the  nose  vertex)  but  for  most  practical  purposes  are  considered  pointed 
bodies.  Blunt-nosed  bodies,  of  course,  violate  the  basic  assumptions  made  in 
the  development  of  the  rule,  i.e.,  slender  bodies  and  hypersonic  flow.  In 
addition,  application  of  the  rule  requires  that  for  similarity  in  pressure  or 
drag,  these  parameters  must  be  a  function  of  K  alone  at  any  point  along  the 
body  surface.  The  pressure  at  the  nose  tip  of  a  blunt  body  is  the  stagnation 
pressure  which  is  solely  a  function  of  Mq  instead  of  K.  Thus  for  truly  blunt 

bodies,  like  the  hemisphere,  values  of  C have  not  been  correlated  on  the  basis 
of  the  hypersonic  similarity  parameter. 


Fig.  2  RANGE  OF  APPLICABILITY  OF  SIMILARITY  LAW  FOR  OGIVES 
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D.  PRESSURE  DISTRIBUTIONS  FOR  NON-LIFTING  "POINTED  BODIES"  OF 
REVOLUTION  -  SHOCK  ATTACHED 

1.  Tangent  Ogives 

The  tangent  ogive  nose  shape  is  used  quite  extensively  in  the  de¬ 
sign  of  high  speed  aircraft  and  guided  missiles.  Tangent  ogiv  s  have  the  ad¬ 
vantage  of  a  larger  nose  vertex  angle  and  higher  volume  than  that  of  an  inscribed 
cone  having  the  same  l^ngth-to-diameter  ratio,  but  the  disadvantage  of  slightly 
increased  drag.  This  disadvantage  is  usually  overshadowed,  however,  particularly 
with  reference  to  those  fixed  length  guided  missiles  which  normally  use  the  nose 
to  house  a  radar  dish  or  scanner  as  far  forward  on  the  missile  as  possible. 

The  analysis  required  for  the  development  of  the  "Present  Method” 
was  initially  formulated  for  the  tangent  ogive  nose  shapes  and  before  proceeding 
with  the  method  description,  a  brief  discussion  of  the  geometric  characteristics 
of  the  tangent  ogive  is  presented. 

The  tangent  ogive  is  a  pointed  (for  R  >  r^)  convex  surface  of 
revolution  generated  by  rotation  of  the  radius  vector  R  to  produce  a  circular 
arc  with  the  tangent  at  the  maximum  radius  (r^)  being  parallel  to  the  axis  of 
symmetry  of  the  body  of  revolution. 


The  following  relations  are  used  to  define  the  geometric  character¬ 
istics  of  the  tangent  ogives: 

a.  Body  Surface  curvature  angle,  G 
sin  9  -  --^-x  ,  G  -  sin-’  ^ 

at  the  nose  vertex,  G„  “  sin 
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b.  Caliber,  N 

N  -  R/d 


7n  terms  of  fineness  ratio, 
1 


N  - 


R 

l 


t 

~T 


(-4-)' 


l./d 


c.  Relation  for  fraction  of  axial  distance  along  nose 
centerline,  x/; 

x//  -  i  .  ♦  r)j 


Once  an  t/d  value  has  been  selected  for  the  tangent  ogive,  the 
above  geometric  relations  can  quickly  provide  the  caliber,  if  desired,  and  the 
variations  in  body  surface  angle,  0. 

This  section  considers  only  pointed  bodies  of  revolution  with 
shock  attached.  The  tangent  ogive  may  or  may  not  have  shock  detachment  depend¬ 
ing  on  the  Mach  number  and  nose  vertex  angle  of  the  ogive.  For  this  reason,  a 
plot  of  free  stream  Mach  number  versus  ogive  semi-vertex  angle,  G„,  has  been 
presented  in  Figure  3  wherein  the  shock  detachment  region  has  been  defined.  An 
additional  plot  of  the  tangent  ogive  geometric  characteristics  (0^,  Vd,  N)  has 
also  been  presented  in  Figure  4.  This  additional  plot  provides  a  rapid  means 
of  determining  the  shock  region  into  which  any  one  tangent  ogive  will  fall. 


The  prescribed  method  of  this  study  for  calculating  pressure 
coefficients,  Cp,  represents  a  combination  of  the  Generalized  Newtonian  thr ory 
and  the  shock-expansion  method.  The  application  of  this  method  begins  with  the 
Generalized  Newtonian  theory,  i.e., 


C  sin2© 
pmax 


(1) 


where  C  -  C  /sin*©v  (2) 

pmax  PN  N 

and  ©  -  local  body  surface  inclination  angle 

©N  ■  nose  semi-vertex  angle 

C  -  pressure  coefficient  at  nose  vertex 
PN 


For  the  special  case  of  tangent  ogives,  a  more  convenient  form  of  the  C 
equation  might  be:  P 


C 

P 


rmax 


sing0 
PN  sina©N 


(1 


sin*©  -  C 


-  x/Oa 


(la) 
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NOAE  SEMI-VERTEX  AN8LC  ,  eN  (OCA.) 

Fig.  3  SHOCK  DETACHMENT  MACH  NUMBER  FOR  BODIES  OF 
REVOLUTION  WITH  CONICAL  TIP  NOSES 
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Fig.  4  TANGENT  OGIVE  GEOMETRIC  CHARACTERISTICS 
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Using  the  appropriate  value  of  nose  semi-vertex  angle, 

N 

enter  the  conical  flow  tables,  such  as  those  given  in  Table  6  of  Appendix  1,  and 

obtain  a  value  of  C  at  a  -  0°  for  a  cone  of  equivalent  semi-vertex  angle. 

PN 

C  ,  a  constant,  can  now  be  determined  and  the  pressure  coefficient  distrlbu- 
pmax 

tlon  along  the  body  surface  can  be  calculated  since  0  is  a  known  quantity  at  any 

x/t  value.  A  plot  of  C  vs.  K  for  tangent  ogives  with  M  as  a  parameter  is 

pmax  ° 

shown  in  Figure  5. 


Fag  .  5  C  FOR  TANGENT  OGIVES 
pmax 

The  local  Mach  number  distribution  along  the  surface  of  the 
body  may  be  obtained  directly  from  Table  18  of  Appendix  3  once  a  value  of  p/pt 
has  been  determined  as  follows: 


where  p  -  p  along  a  body  surface  meridian. 

t  tN 


~T 


I 


Il 
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Now 


PN  '  Po 


(P/Pt)  -  (P0/Pt  ) 
*  N  °  *N 


PN 


N 


so 


?t  CD  +  PQ/qO 
__N  „  PN  °  ° 

qo 


Cp  *  P/Pt 


:PN  +  po/qo 

(p/pt) 
t  N 


-  Po/qo 


(3) 


and 


(Wpt\ 


P/P,  -  <CP  *  p0/,0)  rr-—J i 


N 


(3a) 


where  (p/p.)  is  determined  by  first  using  the  cone  tables  (such  as  Table  7  of 
1  N 

Appendix  1)  to  obtain  a  value  of  Mjj  corresponding  to  a  particular  value  of  ©N 


and  M  ;  (p/p.)  can  then  be  obtained  directly  from  Table  18  of  Appendix  3 

O  L  Jy 

(Prandtl-Meyer  flow  table)  for  the  given  value  of  M„. 

At  some  point  along  the  body,  the  Cp  calculated  by  the 
Generalized  Newtonian  theory  will  offer  poor  agreement  with  the  exact  solutions. 
At  or  near  the  value  of  x/t  where  the  Cp  from  impact  theory  begins  to  deviate 

from  the  exact  solution,  the  shock-expansion  method  is  used  to  extend  the 

0 

solution  by  matching,  at  that  point,  the  surface  pressure  coefficients  and  pre¬ 
serving  the  surface  streamline  total  head. 


The  x/t  point  at  which  the  impact  theory  should  be  terminated 
for  the  case  of  tangent  ogives  was  found  by  correlating  the  Cp  values  obtained 
from  the  Generalized  Newtonian  theory  with  those  from  exact  solutions  and  this 
point  was  found  to  be  a  function  of  the  hypersonic  similarity  parameter,  K.  By 
investigating  a  wide  range  of  K  values,  an  empirical  expression  was  determined 
for  the  x/t  value  at  which  the  two  methods  should  be  matched: 


(x/t)  -  0.15  K  +  0.60  for  K<  2.667 

m 


*An  attempt  was  made  to  match  both  the  pressure  and  pressure  gradient  for  the 
pointed  bodies  but  the  pressure  gradient,  which  is  a  function  of  body  curvature 
angle,  could  not  be  matched  in  the  region  of  small  curvature  changes  considered. 
It  will  be  shown  later,  however,  that  the  pressure  gradient  can  be  matched  on 
the  hemisphere  nose  in  the  region  where  the  two  methods  are  to  be  matched  since 
the  body  curvature  changes  are  large. 


4.T 


( 


*  r' 


This  equation  has  been  plotted  in  Figure  6 


as  os  04  to  ti  u  to  to  to 


Fig.  6  MATCHING  POINT  VALUE  OF  x/t  FOR  TANGENT  OGIVES  AT  ALL 
MACH  NUMBERS  CONSIDERED 

Once  x/l  is  known,  the  value  of  C^  at  which  the  shock-expansion 

■ethod  starts  (C_  )  is  found  fron  equation  la.  In  order  to  continue  the  SEM 

PSEM 

aft  of  the  matching  point,  the  following  procedure  must  be  used: 

a.  calculate  a  value  of  p/pt  at  the  matching  point  by 
using  equation  3a; 

b.  from  the  Prandtl-Meyer  flow  equations  or  tables13 
(Table  18  herein),  determine  the  value  of  v  (flow  de¬ 
flection  angle  from  M  ■  1)  corresponding  to  p/pt  above; 
this  v  is  used  as  vb  in  the  following  step; 

c.  knowing  the  body  surface  angle  at  the  matching  point, 

0  ,  and  at  other  points  downstream,  G,  additional  values 
in 

of  p/Pt  and  M  aft  of  the  matching  point  may  be  obtained 

by  entering  Table  18  with  a  known  value  of  v 

(v-v  +0  -  G  -  v_  +  AG)  and  proceeding  in  steps  of 

m  m  In 

v  ■  v  +  AG  over  the  remainder  of  the  body  surface; 
in 

d.  the  additional  values  of  p/pt  may  now  be  substituted 
into  equation  3  and  Cp  values  aft  of  the  matching  point 
can  be  calculated. 

(To  illustrate  the  application  of  the  "Present  Method",  numerical  examples  have 
been  presented  in  Appendix  1  and  Appendix  2  for  a  tangent  ogive  of  fineness 
ratio  3  at  MQ  ■  2.  The  afterbody  pressure  distribution,  which  is  discussed  in 
Section  F,  has  also  been  included  as  part  of  the  numerical  examples.) 
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Using  the  procedure  outlined  above,  pressure  coefficient  dis¬ 
tributions  have  been  calculated  using  the  present  method  for  five  tangent  ogives 
with  t/d  values  of  1.5,  2,  3,  6,  and  12  and  for  Mach  numbers  of  1.5,  2,  3,  and 
6.  These  results  have  been  compared  with  the  method  of  characteristics,  the 
shock-expansion  method,  and  the  tangent-cone  method  in  Figure  7.  The  agreement 
of  the  present  method  with  exact  solutions  is  exceptional  for  all  values  of  the 
hypersonic  similarity  parameter  K.  (The  only  instance  where  deviation  was 
apparent  was  for  Mq  ■  1.5  at  K  -  0.5  where  the  Cp  values  from  the  present  method 
matched  those  from  the  exact  solutions  only  when  faired.)  Notable  disagreement 
was  exhibited  by  the  tangent-cone  and  shock-expansion  methods  for  Mq  *  1.5  at 
K  -  0.5,  I  ■  2  at  K  ■  1,  and  Mp  -  3  at  K  -  2;  the  remaining  cases  indicated 
markedly  improved  agreement  with  exact  solutions,  i.e.,  the  high  Mach  number  and 
fineness  ratio  cases  approached  the  accuracy  provided  by  the  present  method. 
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Fig.  7  TANGENT  OGIVE  PRESSURE  COEFFICIENT  DISTRIBUTION, 

K  -  0.5,  1.0,  and  2.0 


The  percent  error  in  drag  due  to  the  disagreement  of  the  Cp 
values  from  the  above  three  methods  with  the  Cp  values  obtained  from  exact 
solutions  has  been  presented  in  table  form  below  for  the  lowest  Mach  number  of 
each  K  value.  These  particular  cases  were  selected  because  they  represent  the 
poorest  conditions  for  the  three  methods;  in  spite  of  this,  the  present  method 
yields  near  perfect  agreement  while  the  tangent-cone  and  shock-expansion  methods 
offer  considerable  disagreement  when  each  method  is  compared  with  exact  solutions. 


Table  1.  Percent  Error  in  Drag 


Method 

Jo 

Vd 

_K_ 

%  Error  in  Drag 

Present 

1.5 

3 

0.5 

2 

Shock  Expansion 

tt 

ft 

ft 

-84 

Tangent-Cone 

It 

ft 

f* 

29 

Present 

2.0 

2 

1.0 

0 

Shock  Expansion 

ft 

ft 

ft 

-29 

Tangent-Cone 

If 

f* 

ft 

12 

Present 

3.0 

1.5 

2.0 

0 

Shock  Expansion 

ft 

ft 

ft 

-6 

Tangent-Cone 

ft 

tt 

1* 
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Ordinarily,  a  presentation  of  "percent  error  in  drag"  only  would 
not  be  sufficient  to  define  the  accuracy  to  which  the  method  might  predict  the 
pressure  distribution  over  the  body.  This,  of  course,  refers  to  the  possibility 
of  some  one  method  predicting  higher  pressures  at  the  vertex  and  lower  pressures 
near  the  base  (or  vice  versa)  and  as  a  result,  the  compensating  errors  in  pres¬ 
sure  would  serve  to  lower  the  error  in  drag  which  is  merely  an  integration  of 
the  pressure  distribution  values.  Fortunately,  this  condition  did  not  preva.’l 
for  the  methods  used  as  comparisons;  a  small  compensating  error  was  noted,  how¬ 
ever,  for  the  present  method  at  MQ  ■  1.5  and  t/d  -  3. 

As  mentioned  earlier,  the  results  obtained  using  the  present 
method  were  correlated  on  the  basis  of  the  hypersonic  similarity  parameter,  K. 
Figure  8  illustrates  the  validity  of  this  correlation. 

2.  Secant  Ogives 

Secant  ogive  is  the  name  applied  to  an  ogival  nose  cut  off  in 
length,  t  ,  from  its  related  tangent  ogive  length,  t,  by  passing  a  plane 
normal  to  the  longitudinal  axis  of  the  ogival  shape.  A  family  of  secant  ogives 
related  to  the  tangent  ogive  (see  pages  11  and  12  for  sketch  and  nomenclature) 


Fig.  8  CORRELATION  OF  PRESSURE  COEFFICIENT  DISTRIBUTIONS  ON 
TANGENT  OGIVES  USING  THE  HYPERSONIC  SIMILARITY  PARAMETER  K 


of  caliber  N  and  nose  curvature  angle,  0,  can  be  expressed  in  particular  by  the 
following  equation: 


x//  -  1  -  /Ra  -  (Y..+  ry_  t 

i  s.o. 

These  shapes  (for  equal  fineness  ratios)  offer  high  volume 
approaching  that  provided  by  the  tangent  ogives  with  the  further  advantage  of 
exhibiting  less  drag  than  either  the  related  tangent  ogive  or  inscribed  cone  of 
the  same  fineness  ratio.  Pressure  coefficient  distributions  have  not  been  pre¬ 
sented  for  secant  ogives,  as  such,  since  they  are  merely  forward  segments  of  the 
tangent  ogives  discussed  in  Section  D,  Part  1. 

3.  Triple  Cone-Tangent  Ogive  Combination 

This  composite  nose  design  has  been  considered  for  two  reasons: 
the  nose  shape  is  an  unusual  one  wherein  a  triple  cone  segment  preceeds  a  tangent 
ogive  section  and  secondly,  experimental  data  were  available  for  comparison  pur¬ 
poses.  Experimental  data,  of  course,  provide  an  excellent  basis  for  comparison 
studies  but  often  times  the  reliability  or  availability  of  such  information  is 
quite  limited. 

Figure  9  presents  a  comparison  of  the  C  values  obtained  using 
the  present  method,  the  shock-expansion  method,  and  the  experimental  data.  The 
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Fig.  9  PRESSURE  COEFFICIENT  DISTRIBUTION  ON  A  TRIPLE 
CONE-TANGENT  OGIVE  COMBINATION 


present  method  indicates  favorable  agreement  at  both  Mq  -  2  and  4  whereas  the 
shock-expansion  method  is  quite  poor  at  Mq  ■  2  but  like  the  present  method,  she  ns 
favorable  agreement  at  M  ■  4.  The  failure  of  the  shock-expansion  method  in  pre¬ 
dicting  the  proper  C values  at  Mq  -  2  was  to  be  expected  since  the  value  of  the 
hypersonic  similarity  parameter,  K  -  0.59,  was  far  too  low  for  this  method. 

4.  von  Karman  Minimum  Drag  Nose  Shapes 

The  von  Karman  nose  shapes  are  high  volume-low  drag  nose  shapes 
developed  by  von  Karman  for  a  given  length  and  diameter  for  practical  fineness 
ratios  and  moderate  supersonic  Mach  numbers.  The  nose  shapes  are  defined  by  the 
equation : 

r  -  v/0  -  1/2  sin  20 

where  0  -  cos-'  (1  -  2x/')  . 


These  nose  shapes  have  mathematically  infinite  slopes  at  their 
vertices  yet  for  most  practical  purposes  are  considered  pointed  bodies.  Using 
the  present  method  for  predicting  Cp,  however,  these  noses  require  special  con¬ 
sideration  in  that  they  cannot  be  treated  as  either  blunted  or  pointed  bodies  of 

revolution  and  the  difficulty  arises  in  the  proper  determination  of  C  .  In 

pmax 


an  effort  to  obtain  starting  values  of  C  to  use  with  the  Generalized  Newtonian 

pmax 

theory,  equation  1,  a  systematic  study  was  initiated  whereby  cones  were  fitted  to 

each  von  Karman  nose  shape  tangent  to  the  nose  at  some  \Q/ >  value.  Using  the  cone 

tables  at  a  given  free  sti’eam  Mach  number,  C  was  determined  and  pressure  dis- 

pmax 

tributions  were  calculated  using  the  Generalized  Newtonian  theory  and  compared 
with  the  theoretical  values4.  When  the  present  method  matched  the  theory,  the 


x{,//  and  Mach  number  values  were  noted  and  the  entire  process  repeated  using 

various  other  fineness  ratios  and  Mach  numbers.  Ultimately,  this  study  provided 

the  curve  of  Figure  10  wherein  xq/  versus  K  has  been  plotted.  Accompanying  this 

curve  is  a  plot  of  C  versus  K  as  shown  in  Figure  11.  This  is  a  fictitious 

pmax 


C  value  to  be  used  wit 
pmax 


the  Generalized  Newtonian  equation  (Eq.  1), 


C 

P 


C  sin®0  , 
pmax 


for  calculating  pressures  on  von  Karman  noses. 


^Solutions  may  be  carried  forward  of  *c/'  using  the  G  at  the  desired  point,  x, 

and  the  empirically  determined  value  of  C  ;  the  error  incurred  is  negligible 

pmax 


as  indicated  by  the  results. 
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As  in  the  cases  for  the  previously  discussed  nose  shapes,  the 
Generalized  Newtonian  theory  has  been  extended  using  the  shock-expansion  method. 
The  matching  point  for  the  two  methods  was  determined  analytically  and  has  been 
plotted  as  ( x/l )  versus  K  in  Figure  12.  These  values  differ  somewhat  from  those 

m 

for  the  ogival  shapes  (see  Fig.  6). 


Fig.  10  TANGENCY  POINT  FOR  CONES  FITTED  TO  VON  KAR1IAN  NOSE  SHAPES 


Fig.  12  MATCHING  POINT  VALUE  OF  x/l  FOR  VON  KARMAN  NOSE  SHAPES 
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The  pressure  distributions  calculated  using  the  present  method  were 
compared  with  the  theoretical1 4  values  at  MQ  -  1.5  and  2  (approximately)  for  the  K 
region  in  which  the  theory  is  applicable,  i.e.,  from  K  -  0.55  to  0.922.  Figure  13 
presents  these  comparisons  and  the  agreement  appears  quite  favorable. 
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Fig.  13  PRESSURE  COEFFICIENT  DISTRIBUTION  FOR  VON  HARMAN 
MINIMUM  DRAG  NOSE  SHAPE,  / /d  -  2.10,  2.32,  and  2.72 


The  use  of  the  present  method  has  provided  a  means  of  extending 
the  restricted  K  region  characteristic  of  the  von  Karman  nose  shapes.  The  C 

p 

xc/r,  and  <x/')m  curves  were  extrapolate*,  to  K  values  of  almost  2  and  pressurlT* 
distributions  were  calculated  and  compared  with  experimental  data16  gathered  at 
Mo  "  1,61  *  2*80>  and  4.00,  This  compariSon  is  presented  in  Figure  14  and  the 
agreement  appears  very  good,  indicating  the  useful  range  of  the  K  values  can  now 
be  more  than  doubled  when  applying  the  Solutions  offered  by  the  present  method. 
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5.  Power  Series  Minimum  Drag  Nose  Shapes 


The  minimum  drag  power  series  now .  often  referred  to  as  the  hyper¬ 
sonic  optimum  nose  since  it  closely  approximates  t ho  profile  of  Newton's  hyper¬ 
sonic  optimum  nose,  is  defined  by  the  equation: 


r 


r 


b 


(x/) 


3/4 


These  nose  shapes,  like  the  von  Karman  noses,  have  mathematically  infinite  slopes 
at  the  vertex  and  calculational  procedures  similar  to  those  performed  for  the  von 
Karman  shapes  had  to  be  carried  out  for  the  3/4  power  noses  with  one  notable  ex¬ 
ception,  i.e.,  the  Generalized  Newtonian  theory  predicts  the  pressure  distribution 
very  accurately  along  the  entire  surface  of  the  nose.  Worthy  of  note  in  this  in¬ 
stance  is  the  geometry  of  the  power  series  nose  shapes  which  have  a  finite  slope 
at  the  shoulder  or  base.  This  eliminates  the  need  for  extending  the  theory  using 
the  SEM  since  this  type  of  body  has  no  shallow  angle  section  wherein  the  General¬ 
ized  Newtonian  theory  is  inapplicable.  Thus  only  plots  of  C  and  x  /  were 

pmax 

predetermined  for  these  nose  shapes  and  arc  presented  in  Figures  15  and  16. 
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Pressure  distribution  coefficients  were  calculated  using  the 

present  method  for  three  different  3/4  power  noses  of  fineness  ratio  3,  5,  and 

7,  for  Mach  numbers  from  2  to  7.5,  These  values  were  compared  with  theoretical' 4 

C  values  in  Figure  17  and  once  again,  the  agreement  is  exceptional.  Note  that 
P 

solutions  were  carried  forward  of  the  x  /L  value  using  the  corresponding  C 

C  pmax 

value  in  the  Generalized  Newtonian  theory  and  the  resulting  error  is  negligible. 

6.  Isentropic  Spikes 

Brief  consideration  has  been  given  these  shapes  even  though  they 
are  primarily  used  as  supersonic  diffusers  for  ramjet  engines.  The  function  of 
the  diffuser  is  to  decelerate  the  air  from  its  free  stream  velocity  at  the  ram¬ 
jet  intake  to  a  velocity  at  the  combustor  which  is  compatible  with  the  available 
flame  velocity. 

The  true  isentropic  spike  would  ideally  have  a  total  pressure  re¬ 
covery  factor  of  one  but  the  length  of  such  a  spike  would  make  it  impractical 
for  use.  The  long  needle-like  nose  would,  of  course,  make  it  structurally  un¬ 
sound  and  the  boundary  layer  build-up  along  the  spike  would  destroy  its  effec¬ 
tiveness.  Due  to  these  conditions,  conical  tips  are  usually  attached  forward 
of  the  "isentropic  surfaces"  with  the  end  result  of  some  loss  in  total  pressure 
recovery  attributed  to  the  how  shock-wave.  The  spikes  are  usually  designed 
such  that  the  bow  shock  wave  and  the  Mach  lines  from  the  compression  surfaces 
coalesce  at  a  point  called  the  focal  point.  This  point,  by  design,  usually 
lies  at  or  very  near  the  cowl  lip  of  the  ramjet  engine  depending  on  the  desired 
conditions  at  the  combustor.  Because  of  the  design  mobility  in  focal  point  and 
possible  compromises  in  total  pressure  recovery  for  a  given  condition,  isentropic 
spikes  of  basically  two  designs  are  available:  those  which  yield  a  constant 
total  pressure  recovery  factor  by  merely  varying  the  cone  angle  and  free  stream 
Mach  number  and  secondly,  those  which  have  a  fixed  cone  tip  angle  thereby  causing 
the  total  pressure  recovery  to  decrease  with  increasing  Mach  number.  These  two 
designs'6’17  were  investigated  in  the  Mach  number  range  from  2  to  4  and  pressure 
distributions  were  calculated  using  the  present  method  and  compared  with  exact 
solutions.  The  matching  point  value  of  6  for  the  present  method  was  arrived 


(HEF  14  ) 


at  analytically  and  has  been  presented  in  Figure  18  as  a  function  of  free  stream 
Mach  number.  As  noted  the  matching  point  0  has  been  non-dimensionalized  and 
given  as  ®m/©N  in  order  to  eliminate  the  variable  ©N;  in  effect,  the  matching 
point  ©  has  been  normalized  with  respect  to  ©N  thereby  eliminating  the  need  for 
a  family  of  curves  to  represent  spikes  of  different  conical  tip  angles. 


Fig.  18  RATIO  OF  MATCHING  POINT  ©  TO  NOSE  SEMI-VERTEX  ANGLE  VERSUS 
FREE  STREAM  MACH  NUMBER  FOR  ISENTROPIC  SPIKES 

Figure  19  presents  the  results  of  the  Cp  comparison  between  the 
present  method  and  exact  solutions  for  an  isentropic  spike  with  a  15°  conical 
tip  at  MQ  -  2.5,  3.0,  and  4.0.  The  agreement  appears  good  at  all  Mach  numbers 
and  exceptional  at  MQ  -  4.0.  It  should  be  emphasized  at  this  point,  however, 
that  the  usefulness  of  the  present  method  diminishes  with  increasing  Mach  number 
as  noted  in  Figure  18,  i.e.,  the  value  of  0  /©„  has  reached  1.1  at  M  -  4  in- 
dicating  the  matching  point  lies  very  near  the  conical  tip.  In  essence,  this 
condition  indicates  the  shock-expansion  method  alone  would  be  adequate  (for  most 
engineering  purposes)  in  determining  Cp  at  the  higher  Mach  numbers  for  isentropic 
spikes. 
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E.  PRESSURE  DISTRIBUTIONS  FOR  NON-LIFTING  BLUNT-NOSED  BODIES  OF 
REVOLUTION 

1 .  Hemisphere 

The  hemispherical-type  nose  shape  exhibits  certain  advantages  for 
supersonic  flight  vehicles  that  are  lacking  for  the  generally  accepted  optimum 
pointed-nose  bodies  of  revolution.  Two  possible  features  of  the  hemisphere  type 
nose  shape  are  its  ability  to  readily  accommodate  radome  systems  and  secondly, 
the  reduced  heat-transfer  effects  at  high  speeds  which  are  characteristic  of 
this  type  body.  The  latter  feature,  of  course,  is  of  utmost  importance  for  re¬ 
entry  bodies. 

The  use  of  impact-shock-expansion  theory  in  determining  the  pres¬ 
sure  distributions  for  hemisphere  type  noses  has  been  used  in  Reference  18  for 
very  high  Mach  numbers  (MQ  8s  7.7)  and  shown  to  be  of  great  value.  Reference  18, 
however,  defines  the  matching  point  only  for  MQ  approaching  infinity  whereas 
the  present  analysis  will  consider  variations  in  the  matching  point  which  result 
for  free  stream  Mach  numbers  less  than  approximately  8. 

Solutions  for  the  pressure  distribution  on  the  hemisphere  nose 
are  initiated  using  the  Generalized  Newtonian  theory: 


C  sin9 9 

pmax 


where  C  ■  C  /sin*0M 

pmax  PN  N 

Since  sinaG„  equals  1  at  the  nose  vertex,  C  -  C  .  For  the  blunt-nosed 
N  Pmax  PN 

bodies,  a  detached  shock  appears  at  the  nose  and  the  value  of  C  becomes  the 

PN 

stagnation  value  behind  the  normal  shock.  This  stagnation  value  is  a  function 
of  free  stream  Mach  number  only  (for  y  "  constant)  and  can  be  readily  calculated 
using  the  Rayleigh  formula: 


For  convenience,  this  equation  has  been  plotted  in  Figure  20  and  as  illustrated 
reaches  a  limiting  value  of  1.84. 

At  some  point  along  the  surface  of  the  hemisphere,  the  Generalized 
Newtonian  theory  begins  to  show  marked  disagreement  with  comparison  data  and  it 
is  at  this  point  that  the  shock-expansion  method  (which  utilizes  the  Prandtl-Meyer 
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Fig.  20  C  FOR  HEMISPHERE  OR  HEMISPHERE-CONE 
*max 

relations)  is  used  to  extend  the  Impact  theory.  This  is  accomplished  by  matching 
the  Cp  values  for  the  two  methods  at  the  point  of  deviation;  this  point  is  de¬ 
termined  by  equating  the  pressure  and  pressure  gradient  given  by  the  Generalized 
Newtonian  theory  and  the  Prandtl-Meyer  relations  as  follows; 

From  Generalized  Newtonian  theory, 

C  -  C  sln'o  -  — - - — 

p  Pmax  (y/2 )po*0* 

Solving  for  p 

p  -  pQ  +(r/2)p  M  •  C  sina0 

^max 

ai  *(>'/2>>>o"o‘  Cp  =»n  28 

*max 

(y/ 2)M  *  C  sin  20 

and  1  dp  _  _ pmax _ 

P  **  l+(y/2)M  8  C  sin*© 

pmax 


3*'  JT 
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From  the  Prandtl-Meyer  relations, 


J_  dp.  rM* 

p  35  y5*-rr 


By  equating  the  above  product  of  -i—  from  the  two  methods,  the  matching  point 
O  can  be  determined.  It  is  apparent  that  the  value  of  0  at  the  matching  point  is 
a  function  of  MQ  and  a  plot  of  Mq  versus  matching  point  0  has  been  presented  in 
Figure  21.  (Note  that  for  very  low  Mach  numbers,  there  is  no  point  at  which  both 
the  pressure  and  pressure  gradient  are  the  same  for  the  two  methods  and  the  curve 
has  been  extrapolated  in  this  region.)  For  sufficiently  high  Mach  numbers 
(M  >21),  the  matching  point  0  reaches  a  limiting  value  of  approximately  35.4°. 
This  fact  is  borne  out  by  Vaglio-Laurin  and  Trella18  as  they  indicate  the  match¬ 
ing  point  Q  for  MQ  approaching  infinity  can  be  obtained  by  the  following  ex¬ 
pression: 


Once  the  matching  point  0  has  been  determined,  the  remainder  of  the  solution  is 
quite  simple  since  p/pt  can  be  calculated  using: 


Cp  "  P/Pt 


CPN  +  Po/qc 
(p/Pt> 

L  VT 


*  Po/qo 


Knowing  p/pt*  the  Prandtl-Meyer  relations  or  tables  provide  the  local  Mach 
number  and  flow  deflection  angles  (from  M  -  1)  so  that  one  may  proceed  along  the 
surface  of  the  body  according  to  v  +  AO  (see  page  16) . 

Pressure  coefficient  distributions  were  calculated  for  the  hemi¬ 
sphere  nose  using  the  present  analysis  and  compared  with  experimental  data1  9  at 
MQ  ■  1.82,  2.81,  3,74,  and  4.76  and  a  theoretical  treatment  using  numerical 

analysis10  at  MQ  *•  7.7.  The  agreement  of  the  present  analysis  with  the  com¬ 
parison  data  is  quite  favorable  as  illustrated  in  Figure  22. 

2.  Hemisphere -Cone  Combination 

The  hemisphere-cone  type  nose  has  found  extensive  use  among  high 
speed  vehicles  for  much  the  same  reasons  as  the  hemisphere  nose.  The  hemisphere- 
cone  nose,  although  not  exhibiting  the  high  volume  of  a  hemisphere,  does  have 
the  advantage  of  reducing  the  heat-transfer  effects  at  high  speeds  and  in  addition 
provides  somewhat  more  stability  than  that  of  the  hemisphere.  Recent  studies  by 
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Fig.  21  MATCHING  POINT  VALUE  OF  0  FOR  HEMISPHERE  OR 

HEMISPHERE-CONE 


Fig.  22  HEMISPHERE  PRESSURE  COEFFICIENT  DISTRIBUTION 
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Perkins,  Jorgensen  and  Sommer’0  have  also  indicated  that  the  drag  of  a  nose 
shape  consisting  of  a  hemispherical  surface  faired  into  an  expanding  conical 
surface  can  be  less  than  that  of  a  sharp  cone  of  the  same  fineness  ratio  at  all 
supersonic  Mach  numbers. 

Pressure  distributions  have  not  been  calculated  for  the  hemi¬ 
sphere-cone  nose  shapes  since  the  procedure  would  be  similar  to  that  of  the 
hemisphere.  Moreover,  the  present  analysis  will  not  predict  the  over-expansion 
of  the  flow  at  the  hemisphere-cone  juncture  which  occurs  for  Mach  numbers  less 
than  about  3.  Above  Mq  -  3,  however,  the  present  analysis  should  predict  the 
pressure  coefficients  accurately  enough  for  engineering  needs  since  the  expanding 
flow  over  the  hemispherical  tip  does  not  reach  a  lower  pressure  than  the  exact 
value  for  a  cone  having  the  same  slope  as  the  conical  afterbody  and,  in  addition, 
this  pressure  remains  essentially  constant  over  the  entire  conical  afterbody  as 
would  be  predicted  by  impact-shock  expansion  theory. 

Reference  20  presents  experimental  data  which  illustrate  the  flow 
over-expansion  phenomenon  associated  with  the  hemisphere-cone  type  nose  shapes 
for  Mach  numbers  less  than  3.  The  extent  of  the  over-expansion  as  pointed  out 
in  this  reference  is  dependent  primarily  on  the  slope  of  the  conical  afterbody 
surface;  as  the  slope  decreases,  the  over-expansion  of  the  flow  appears  more 
prominent.  At  all  Mach  numbers,  however,  the  pressure  on  the  conical  afterbody 
returns  almost  to  the  exact  value  of  a  cone  with  a  slope  corresponding  to  that 
of  the  conical  afterbody.  With  this  in  mind,  a  first  approximation  for  pre¬ 
dicting  the  pressure  coefficients  of  hemisphere-cone  noses  by  the  present  anal¬ 
ysis  for  Mq  less  than  3  is  quite  possible.  It  has  been  shown  that  the  present 
analysis  would  give  accurate  results  over  the  hemisphere  portion  of  the  nose  and 
since  experimental  data  have  indicated  the  conical  afterbody  pressure  approaches 
that  of  the  equivalent  cone,  it  would  be  within  the  acceptable  accuracy  limits 
of  most  engineering  methods  to  merely  fair  the  present  analysis  solution  into 
the  equivalent  cone  solution. 

F.  PRESSURE  DISTRIBUTIONS  ON  NON-LIFTING  CYLINDRICAL  AFTERBODIES 

Fenter's13  method  has  been  used  to  predict  the  pressure  distributions 
over  cylindrical  afterbodies.  This  method,  which  is  a  simplification  of  the 
second-order  shock  expansion  method,  yields  the  following  relation  for  pressure 
distribution : 

s 

Pn/po  "  pcn/po  +  (pn'/po  '  Pcn^o’  e  ^ 

where  s  -  segment  length 
l  —  nose  length 
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K  ■  hypersonic  similarity  parameter,  V^d// 
d  ■  diameter  of  body 

with  subscripts  and  superscripts 

n  ■  specific  segment  of  equivalent  tangent  body 

cn  “  condition  c;:  the  equivalent  tangent  cone  of  a  specific  segment 

'  -  condition  at  the  most  forward  point  of  an  equivalent  tangent 
body  segment  (s  -  0) 


y 


The  values  of  Pn/PQ  and  Pcn/PQ  are  the  actual  pressures  on  the  ogive 
(tangency  points)  whereas  the  primed  values  pn'/po  represent  conditions  at  the 
most  forward  point  of  any  segment  and  are  obtained  using  the  Prandtl-Meyer  re¬ 
lations  or  tables  for  expansions  through  ©R  -  ®n+j.  etc. 

Since  the  Fenter  method  is  being  used  on  the  cylindrical  afterbody 
alone,  the  relation  for  pressure  distribution  reduces  to: 


(4) 


where: 

Ax  - 

segment  length 

l  - 

nose  length 

and 

Pen  _ 

1  since  the  equivalent  tangent  cone  angle 

Pq 

and  C  -0 

pcn 

0° 


Using  equation  4,  the  pressure  coefficient  becomes: 


(5a) 


or 


Ax// 


(5b) 


where  Ax//  -  non-dimensional  distance  from  base  of  nose,  where  x/t  -  1,  to 
the  desired  point  on  the  afterbody, 

G.  APPLICATION  OF  PRESENT  METHOD  TO  LIFTING  POINTED  BODIES  OF  REVOLUTION 

Reliable  theory  amenable  to  rapid  hand  calculatlonal  procedures  for  de¬ 
scribing  the  pressure  distributions  along  lifting  bodies  of  revolution  is  prac¬ 
tically  non-existent,  especially  for  higher  angles  of  attack.  The  pressure 
distributions  provided  by  the  present  method  for  angles  of  attack  up  to  10  de¬ 
grees  yield  results  which  are  comparable  and  at  times  better  than  various  other 
available  methods.  Theoretical  treatment  of  the  problem  at  higher  angles  of 
attack  is  hampered  by  the  effects  of  cross-flow  separation  on  the  leeward  side 
of  the  body  and  compressibility  effects  when  at  high  cross-flow  Mach  numbers. 

In  the  following  two  sub-sections,  the  application  of  the  present 
method  to  bodies  of  revolution  at  low  and  moderate  angles  of  attack  is  described 
and  numerical  examples  are  presented  in  the  appendices. 

1.  Bodies  of  Revolution  at  Small  Angles  of  Attack 

The  present  method  for  predicting  pressure  distributions  on  non-lifting 

bodies  will  also  be  used  on  the  lifting  bodies  of  revolution.  The  combination 

of  Generalized  Newtonian  theory  and  the  shock-expansion  method  will  be  used  to 

define  C  for  the  nose  section  and  the  method  of  Fenter1*  will  be  used  to  con- 
P 

tinue  the  solution  over  the  cylindrical  afterbody.  Beginning  with  the  Generalized 
Newtonian  theory,  Equation  1  is  modified  to  read: 

C  -  C  sina6  (6) 

p  pmax 


where  C  -  C  /sin*  6 

pmax  pN 


(7) 


sin  5  -  sin  O  cos  a  -  cos  *  cos  0  sin  a  (8) 

and  0  -  body  surface  inclination  angle 

<|i  -  radial  meridian  along  body  surface 
a  -  angle  of  attack 

E>  -  angle  between  free  stream  and  tangent  to 
the  body  surface 


O 


For  a  -  0°,  sin  6  **  sin  0. 


The  value  for  C  on  the  desired  meridian,  *,  may  be  determined 
Pjnax 

either  experimentally  or  theoretically.  The  NASA*1 ,aa  zero  and  small  angle 

M 

of  attack  cone  tables  can  be  used  for  the  starting  values  on  a  cone  fitted  to 
the  nose  vertex  if  experimental  data  are  lacking;  the  procedure  is  as  follows: 
The  general  expression  for  pressure  coefficient  is: 


JL 

Po 


JL 

P. 


Pi 


where  p,  refers  to  conditions  at  zero  angle  of  attack  and  the  quantity  p, /pQ 
may  be  obtained  directly  from  the  zero  angle  of  attack  cone  tables**  .  The  ratio 
of  p/p,  which  is  the  ratio  of  static  pressure  at  angle  of  attack  to  static 
pressure  at  zero  angle  of  attack  can  be  calculated  using  the  theory  of  Stone*3 
wherein  the  velocity,  pressure,  and  density  are  expanded  in  the  following  series 
and  higher  order  terms  in  a  are  neglected: 

M*  -  vC  -  cos  t 
p  -  Pi  +  aps  cos  t 

p  ■  p,  +  opt  cos  * 

where  the  flow  quantities  M,*,  p,  ,  and  o,  refer  to  conditions  at  zero  angle  of 

attack  and  l£,  p2 ,  and  os  are  the  flow  quantities  related  to  the  effect  of  angle 

of  attack.  Reference  23  provides  solutions  for  the  above  equations  yielding, 


*The  necessary  tables  required  in  the  calculational  procedures  of  this  sub¬ 
section  have  been  reproduced  from  the  references  and  included  in  Appendix  1. 


among  others,  the  following  expression: 


References  21  and  22  have  provided  tabulated  values  of  II *,  It, 

and  S-  /R> for  use  in  determining  p/p;  and  subsequently,  C^.  Since  conditions  at 

the  nose  vertex  were  used  to  calculate  the  pressure  coefficient,  this  value  of 

C  becomes  C  from  whence  C  ,  Equation  (7),  may  be  determined.  Starting 
p  PN  pmax 

values  of  C  must  be  calculated  for  each  value  of  *  at  the  desired  angle  of 
pmax 

attack.  Once  this  is  done,  the  Generalized  Newtonian  theory,  Equation  (6),  is 
used  to  calculate  the  values  of  along  each  meridian  up  to  the  point  where  the 
shock  expansion  method  is  used  to  extend  the  impact  theory;  for  the  cylindrical 
afterbody,  Equation  (5a)  or  (5b)  is  used  for  the  pressure  coefficient  distribut¬ 
ion. 

The  1  >cal  Mach  number  at  the  nose  vertex  (where  M  -  V  may  be 
determined  from  l  \0  -elation84 : 


M8 


(y+1)  -  (y-1)  M*s 


On  the  leeward  side  of  the  body,  certain  portions  of  the  surface 
will  lie  in  the  "aerodynamic  shadow"  and  the  Newtonian  theory  cannot  predict  the 
pressures  in  this  region.  This  shadow  area  on  the  body,  shown  pictor tally  here. 


can  be  determined  quite  easily  along  any  leeward  meridian  by  setting  sin  h  -  0 
and  solving  for  0.  For  any  meridian  in  general,  the  0  at  which  the  shadowed 
area  starts  is: 

0  ■  t<an-1  (tan  a  cos  <■) 
for  *  -  0°  ©  ■  a 


Thus  the  shadow  area  begins  along  the  t  -  0°  meridian  at  the  point  where  the 
body  surface  angle  is  equal  to  the  angle  of  attack. 

On  all  meridians,  the  matching  point  (where  the  shock-expansion 
method  is  used  to  extend  the  impact  theory)  is  obtained  directly  from  plots 
presented  in  prior  sections  ^f  this  report,  e.g.,  Fig.  6,  page  16.  When  the 
matching  point  lies  within  the  shadowed  area,  the  ^nock-expansion  method  should 
be  started  at  the  x/'  value  where  sin  6  ■  0. 

A  numerical  example  is  presented  in  Appendix  1  wherein  the  pressure 
coefficient  distribution  has  been  calculated  for  a  tangent  ogive-cylindrical 
afterbody  combination  at  Mq  -  2.0  and  t/d  -  3  at  an  angle  of  attack  of  5  degrees. 
Tables  are  included  which  provide  the  necessary  parameters  for  determining  the 
pressure  and  local  Mach  number  distribution.  For  the  numerical  example,  start¬ 
ing  values  of  C  were  determined  for  seven  meridians  beginning  at  if  -  0°  and 
PN 

proceeding  in  30  degree  increments  up  to  *  -  180°.  The  pressure  distribution 
along  the  t  ■  180°  meridian  is  the  only  one  for  which  the  Cp  calculation  is  shown 
since  the  procedure  is  similar  for  each  meridian.  The  calculation  of  the  cylin¬ 
drical  afterbody  pressure  distribution  using  Equation  (5b)  has  also  been  in¬ 
cluded  for  the  -  180°  meridian. 

The  pressure  distributions  calculated  by  the  present  method  for  the 
numerical  example  of  Appendix  1  are  presented  in  Figure  23  where  experimental 


Fig.  23  PRESSURE  COEFFICIENT  DISTRIBUTION  FOR  A  TANGENT 
OGIVE-CYLINDRICAL  AFTERBODY  COMBINATION  AT  a  ■  5", 
i  -  0°,  30°,  60°,  90°,  120°,  150°,  and  180° 


values9*  have  been  included  for  comparison  purposes  with  the  subsequent  indi¬ 
cation  of  good  agreement. 

2.  Bodies  of  Revolution  at  Large  Angles  of  Attack 

The  calculation  of  pressure  distributions  for  angles  of  attack 
greater  than  about  5  degrees  requires  somewhat  more  time  since  the  starting 

values  of  C  must  be  determined  from  the  MIT8®’*7’98  cone  tables  (or  comparable 
PN 

information)  unless  experimental  data  are  available.  Once  the  starting  values 
have  been  obtained,  the  procedure  is  similar  to  that  used  for  the  small  angles 
of  attack. 

Since  the  MIT  cone  tables  are  being  used  to  obtain  starting  values 
at  the  nose  vertex,  the  symbols  and  nouenclature  adopted  by  this  widely  used 
reference  will  not  be  altered  to  conform  to  the  nomenclature  of  this  report. 

Any  attempt  to  redefine  or  re-reference  the  parameters  would  inevitably  tend  to 
compound  the  existing  complexity  of  the  equations.  The  symbols  and  nomenclature 
from  the  MIT  cone  tables  which  are  used  in  the  present  study  are  defined  on  page 

4  of  this  report.  In  addition,  the  parameters  required  for  calculating  C  have 

PN 

been  picked  from  the  MIT  cone  tables  and  are  tabulated  in  Appendix  2  for  ready 
use.  If  any  intermediate  values  of  the  parameters  are  required,  they  can  be 
obtained  much  more  readily  when  nomenclature  consistent  with  the  reference  source 
is  used. 

The  theory  of  Stone93 is  used  once  again  to  determine  the  flow 
parameters  which  have  been  expanded  in  the  following  series;  for  large  angles 
of  attack,  the  higher  order  terms  in  a  cannot  be  neglected,  i.e., 

P/p"  -  1  +  a  A,  cos  *  +  a8  (As  +  A.-,  cos  2#) 

o/F  -  1  +  a  B,  cos  ♦  +  a8  (Bs  +  B3  cos  2*) 

where  p/p  and  p/F  are  the  pressure  and  density  on  the  cone  surface  at  angle  of 
attack  divided  by  the  corresponding  values  at  zero  angle  of  attack,  A,  and  B, 
specify  the  first  order  effects  of  a,  and  Aa  ,  A3  ,  B-,  ,  and  B3  represent  the  second 
order  effects  of  a.  Proceeding  further: 

A.  -  -Vp 

As  -  pQ/p  +  ~  +  ZjT  Cot  °s 

As-  p- /p  +  - - —  cot  e 

Z  0s  2p  s 


B,  -  -l/o 


Bg  m  p  / p  + 


1  u8  .  % 


_ _  + - cot  0 

5a  2o  8 


Bj  ■  os /o  + 


cot  0. 


By  entering  the  cone  tables  at  the  desired  free  stream  Mach  number 
(corresponds  to  M  with  no  subscript  in  tables)  and  the  nose  semi-vertex  angle 
©N  equal  to  ©g,  the  following  quantities  may  be  obtained: 

*  u/c  and  a*/c* 

VP  and  */p 

PQ/P.  Pj/p.  pq/p  and  pa /c 

With  these  quantities,  A,  ,  As  ,  A3  ,  B,  ,  B3  ,  and  B3  may  be  evaluated 
and  subsequently  p/p  can  be  obtained.  This  value  of  p/p  must  now  be  referred 
to  the  free  stream  in  order  to  calculate  the  surface  pressure  coefficient,  i.e., 

P  .  P  _p_  (subscript  1  denotes  free 
Pi  p  p,  stream  in  tables) 

-  p 

where  and  the  quantities  p/p  and  p  /p,  are  taken  directly  from  the 

tables  herein.  The  quantity  p/p  will  be  listed  as  p_/p„  (p_  is  the  cone  surface 

W  B  W  B 

static  pressure  and  pw  is  the  static  pressure  immediately  behind  the  shock  wave). 

Knowing  p/p,  ,  the  pressure  coefficient  at  any  meridian,  *,  on  the 
cone  surface  can  be  calculated  using: 


cp-p'/'"(lr-  *) 


Since  the  value  of  Cp  for  a  given  *  is  determined  at  the  nose  vertex,  this 
represents  the  starting  value  of  C 

PN 

To  compute  the  local  Mach  number  at  the  nose  vertex,  we  begin  by 
expressing  the  entropy  in  powers  of  a  similar  to  that  for  pressure  and  density: 

S-S-aS,  +a*  (Ss  +  S9 ) 


These  values  appear  as  u  and  a8  in  Reference  26. 


I 
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where  S:  -  Cy  [a,  -  y  Bi] 

S,  -  %  [y  B,a/4  -  A*/4  +  Ae  -  y  Bs] 

S,  -  cy  [y  B*  /4  -  Af/4  +  A3  -  y  B3J 
Now  from  the  equation  of  state, 

-  (S-S) 

o/p  -  (p/p)1/y  e  ycv 

p/p  has  already  been  determined  and  (S- SO/yCy  may  be  calculated  using  the  afore¬ 
mentioned  entropy  distribution  equations  wherein  the  quantities  A,  ,  As ,  A3 ,  B,  , 

B2  ,  and  B3  have  also  been  determined  (p/p  can  also  be  calculated  using  the 
equation  on  page  40),  Finally,  from  Bernoulli's  equation, 

V*  +  “  c*  wh^re  c  -  limit  speed. 

pTy-IT 

Dividing  by  a*  and  simplifying, 

-  £l  _  2  _  ca  a*  2 

"  a’  ’  r-i  "  i»  ,i  ‘  Y- 1 

where:  —  -  E^E  and  ca/aa  is  obtained  from  Table  17  herein. 

a*  c/p 

A  numerical  example  has  been  provided  in  Appendix  2  which  illus¬ 
trates  the  manner  in  which  Cp  was  obtained  for  a  tangent  ogive-cylindrical  after¬ 
body  combination  at  Mq  ■  2.0,  Vd  -  3  and  a  -  10°.  As  was  done  for  the  numerical 

example  at  5  degrees,  starting  values  of  C  were  calculated  and  are  shown  for 

PN 

v  ■  0°  to  180°  in  30°  increments.  Only  the  calculation  for  Cp  values  along  the 
w  -  180*  meridian  (including  the  pressure  distribution  on  the  cylindrical  after¬ 
body)  are  shown  in  the  numerical  example. 

Figure  24  presents  the  pressure  distribution  coefficients  calculated 
by  the  present  method  and  experimental  values*6  which  were  used  for  comparison 
purposes.  The  agreement  at  a  -  10°  is  not.  as  good  as  that  for  the  a  -  5°  case 
especially  on  the  leeward  meridian  but  this  was  more  or  less  to  be  expected  since 
flow  separation  may  now  be  present.  The  more  windward  meridians  indicate  good 
agreement  with  the  experimental  values. 


Fig.  24  PRESSURE  COEFFICIENT  DISTRIBUTION  FOR  A  TANGENT 
OGIVE-CYLINDRICAL  AFTERBODY  COMBINATION  AT 


APPENOtX  I 


CALCULATION  OF  PRESSURE  COEFFICIENT  DISTRIBUTION  FOR 
A  TANGENT  OGIVE- CYLINDRICAL  AFTERBODY  COMBINATION 
AT  AN  ANGLE  OF  ATTACK  OF  S  OEGREES 


A.  CONDITIONS 

FREE  STREAM  MAON  NUMBER  —  M.«  t 

NOSE  FINENESS  RATIO  -  i/d*S 

NOSE  SEMI- VERTEX  ANSLE  -  ©N  •ISSCt* 

HYPERSONIC  SIMILARITY  PARAMETER  -  K».S7 

RATIO  OF  SPECIFIC  HEATS  —  Y*  1.406 

a  PRESSURE  COEFFICIENT  AND  MACH  NUMBER  STARTING  VALUES ,  X//  >0 

i£/B#MP/p#-|)*<p#/\)[wi{ia»/i»  )-iJ  eo.  ai-i 


y 

P/P,  *  Z*  i  j  * It  7  .‘*<8, /*)«»*  ».  Al-t 


M*«  «•*-  «mJ  cost, 


a  IN  RADIANS 


EO.  Al-S 


m  r  «■** 

[|y»0-(r-c)M* 


1 


EO.  AI-4 


WHERE 

•  I.SSSS 

FROM  APPENDIX  1  , 

TABLE  B, 

PAOE 

so 

•c 

•  L424I 

N 

M 

M 

TABLE  s. 

• 

_il_. 

"I 

•  -.7E8S 

■ 

• 

M 

TABLE  K>, 

• 

62  _ 

St/R».OSSi 

M 

M 

M 

TABLE  II, 

m 

83  _ 

*/q#  -3971.  ~ 


NUMERICAL  SOLUTIONS  FOR  EQUATIONS  AI-I, Al-E,  Al -S.AN0  AI-4  ARE  SHOWN  IN 
OOLUMNS  S.S  ,  4,  AN0  S  RESPECTIVELY  OF  TABLE  I,  RASE 


TABLE  2.  PRESSURE  COEFFICIENT  AND  MACH  NUMBER  STARTING  VALUES ,  X//«0,  a* 5 


APPENDIX  I 


D.  PRESSURE  COEFFICIENT  AND  LOCAL  MACH  NUMBER  DISTRIBUTIONS 
ALONG  f  ■  ISO*  MERIDIAN  (SEE  FISURE  IS) 


CONDITIONS  AT  NOSE  VERTEX 


Mm 

■  I.4SS 

— 

OOLUMN 

s. 

TABLE  t 

•.ET40 

- 

• 

10, 

M 

C'N 

■  .4117 

- 

M 

1. 

M 

CPL..„ 

•  t.soss 

— 

EO.  7 

l-C,  VALUES  (COLUMN  I  .TABLE  4)  PROM  THE  NOSE  VERTEX  UR  TO  THE  MATCH INO 
POINT,  (X//^,  ARE  OBTAINED  UMN4  EQ.S.  (THE  MATCH  I  N«  POINT  SOLUTION  IS 
OISCUSSEO  ON  RASE  IS). 

tCp  VALUES  (COLUMN  ♦, TABLE  4)  PROM  (X//»m  TO  */!•  I  ARE  CALCULATED  USINS 
THE  RROCEOURE  OUTUNED  ON  RASES  IS  AND  IS  .FOR  THIS  PARTICULAR  EXAMPLE, 
Cp  »t  BOSS  (P/Rt»-  .8871  (EOS) 

BALL  LOCAL  MACH  NUMBER  VALUES  (OOLUMN  S, TABLE  4)  ARE  OBTAINED  FROM  TABLE 
M,  ARRENOIX  S,  USINS  THE  APPROPRIATE  VALUE  OF  P/Pf  (COLUMN  S,  TABLE  4) 
CALCULATED  FROM  EO.  SO. 


TABLE  4.  PRESSURE  COEFFICIENT  DISTRIBUTION  ANO  MACH  NUMBER 
ALONG  ||f» ISO*  MERIDIAN,  a ■  5* 


GN  ■  OENERALIZEO  NEWTONIAN 
SEM  *  SHOCK  EXPANSION  METHOD 

G  •  BOOT  SURFACE  ANSLE 

.7m  •  MATCH  INS  POINT  OF  SENERALIZEO  NEWTONIAN  ANO  SHOCK  EXPANSION  METHOD 

M  •  M  OORRESPONDINS  TO  P/Pf  FROM  PRANOTL  -  METER  TABLES 

V  •  PRANOTL -MEYER  FLOW  DEFLECTION  ANSLE  OOARESPONOtNS  TO  M 
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APPENDIX  I 

E.  CYLINDRICAL  AFTERBODY  PRESSURE  COEFFICIENT  DISTRIBUTION 
ALONG  1^*180*  MERIDIAN 


U8H 


C  i  C 

P  >  (X fJL»  I ) 


AXAt 


VALUES  OF  Cp  ON  THE  CY'  «DRICAL  AFTERBODY  ARE 
CALCULATED  AS  SHOWN  IN  TABLE  S. 


CpOt/AM)*  _  os56>  •'A0**  COLUMN  4, TABLE  *. 


TABLE  S.  PRESSURE  COEFFICIENT  DISTRIBUTION  ON  CYLINDRICAL 
AFTERBODY,  0« 5* 


1 

8 

3 

4 

x/l 

AX/A 

•  K 

Cp*-.033«  © 

1.0 

0 

1.0000 

-.0338 

l.l 

.1 

.8807 

-.0288 

l.t 

.« 

.7408 

-.0248 

1.3 

.3 

.8878 

-.0214 

1.4 

.4 

.8488 

-.0184 

I.S 

.8 

•47E4 

-.0188 
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*  APPENDIX  t 

CALCULATION  OF  PRESSURE  COEFFICIENT  DISTM1UTION  FOR 
A  TANGENT  OGIVE- CYLINDRICAL  AFTERBODY  COMBINATICN 
AT  AN  ANGLE  OF  ATTACK  OF  10  DEGREES 


A.  CONDITIONS 

FREE  STREAM  MAOH  NUMBER 
NO«e  FINENESS  RATIO 
NOSE  SEMI- VERTEX  ANSLE 
HVFERSONIO  SIMILARITY  PARAMETER 
RATIO  OF  SFECIFIO  MEATS 
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NUMERICAL  SOLUTIONS  FOR  EQUATIONS  At-I  ,  AND  AS-t  ARE  SHOWN  IN 
COLUMNS  S  ANO  7  RESPECTIVELY  OF  TASLE  It,  FARE  89. 
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THE  NOMENCLATURE  USEO  IN  AFFENOIX  X  CORRESPONDS  TO  THAT  USED  IN 
THE  MIT  CONE  TABLES  OF  REFERENCES  SS.S7, ANDES. 
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TABLE  12.  PRESSURE  COEFFICIENT  STARTING  VALUES,  XAf  »0, « » 10 
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VALUES  OF  WT>I  ARE  OS  TAMED  FROM  TABLE  IS  USWtS 
IN  COLUMN  S.  TABLE  IS. 


APPENOIX  2 


E.  PRESSURE  COEFFICIENT  AN 0  LOCAL  MACH  NUMBER  DISTRIBUTIONS 
ALONG  f » ISO*  MERIDIAN  (SEE  FIGURE  >4) 

THE  MME  PROCEDURE  OUTLINED  FOR  THE  EXAMPLE  OF  APPENOIX  I  HA*  SEEN 
USED  IN  THE  FOLLOW  I  NO  EXAMPLE  TD  OBTAIN  THE  PREBBURE  ANO  LOCAL 
MACH  NUMBER  DISTRIBUTION*  ON  THE  NOSE  ANO  CYLINDRICAL  AFTER  BOOT  i  THE 
RESULT*  HAVE  BEEN  SUMMARIZED  IN  TABLES  IS  ANO  IB. 

CONDITIONS  AT  NOSE  VERTEX 

Mu  *  LSS4  -  COLUMN  B,  TABLE  IS 

(P/Pt)M  *.»30B  -  “  S, 
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CL  •  t.tSSB  EO.  7 
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TABLE  IS.  PRESSURE  COEFFICIENT  DISTRIBUTION  ANO  MACH  NUMBER 
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V  •  PRANOTL- MEYER  FLOW  DEFLECTION  ANSLE  CORRESPOND  MO  TO  M 
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APPENDIX  3 

TABLE  18.  PRANDTL- MEYER  FLOW  PARAMETERS 


M 

P/Pt 

V 

M 

P/Pt 

V 

1.00 

0.5283 

0 

6.00 

.6334  -3 

84.955 

1.10 

.4684 

1.336 

6.10 

.5721  -3 

85,635 

1.20 

.4124 

3.558 

6.20 

.5173  -3 

86.295 

1.30 

.3609 

6.170 

6.30 

.4684  -3 

86.937 

1.40 

.3142 

8.987 

6.40 

.4247  -3 

87,561 

1.50 

.2724 

11.905 

6.50 

.3855  -3 

88.169 

1.60 

.2353 

14.861 

6.60 

.3503  -3 

88.760 

1.70 

.2026 

17.810 

6.70 

.3187  -3 

89.335 

1.80 

.1740 

20.725 

6,80 

.2902  -3 

89.895 

1.90 

.1492 

23 . 586 

6.90 

.2646  -3 

90.441 

2.00 

.1278 

26.380 

7.00 

.2416  -3 

90.973 

2.10 

.1094 

29.097 

7.10 

.2207  -3 

91.492 

2.20 

#.9352  -1 

31.732 

7.20 

.2019  -3 

91.997 

2.30 

.7997  -1 

34.283 

7.30 

.1848  -3 

92.490 

2.40 

.6840  -1 

36.746 

7.40 

.1694  -3 

92.971 

2.50 

.5853  -1 

39.124 

7.50 

.1554  -3 

93.440 

2.60 

.5012  -1 

41.415 

7.60 

.1427  -3 

93.898 

2.70 

.4295  -1 

43.621 

7.70 

.1312  -3 

94.345 

2.80 

.3685  -1 

45.746 

7.80 

.1207  -3 

94.782 

2.90 

.3165  -1 

47.790 

7.90 

.1111  -3 

95.208 

3.00 

.2722  -1 

49.757 

8.00 

.1024  -3 

95.625 

3.10 

.2345  -1 

51.650 

8.10 

.9449  -4 

96.032 

3.20 

.2023  -1 

53.470 

8.20 

.8723  -4 

96.430 

3.30 

.1748  -1 

55.222 

8.30 

.8060  -4 

96.820 

3.40 

.1512  -1 

56.907 

8.40 

.7454  -4 

97.200 

3.50 

.1311  -1 

58.530 

8.50 

.6898  -4 

97.573 

3.60 

.1138  -1 

60.091 

8.60 

.6390  -4 

97.937 

3.70 

.9903  -2. 

61.595 

8.70 

.5923  -4 

98.293 

3.80 

.8629  2 

63.044 

8.80 

.5494  -4 

98.642 

3.90 

.7532  -2 

64.440 

8.90 

.5101  -4 

98.984 

4.00 

.6586  -2 

65.785 

9.00 

.4739  -4 

99.319 

4.10 

.5769  -2 

67.082 

9.10 

.4405  -4 

99.646 

4.20 

.5062  -2 

68.333 

9.20 

.4099  -4 

99.967 

4.30 

.4449  -2 

69.541 

9.30 

.3816  -4 

100.282 

4.40 

.3918  -2 

70.706 

9.40 

.3555  -4 

100.590 

4.50 

.3455  -2 

71.832 

9.50 

.3314  -4 

100.892 

4.60 

.3053  -2 

72.919 

9.60 

.3092  -4 

101,188 

4.70 

.2701  -2 

73.970 

9.70 

.2886  -4 

101.479 

4.80 

.2394  -2 

74.986 

9.80 

.2696  -4 

101.763 

4.90 

.2126  -2 

75.969 

9.90 

.2520  -4 

102.043 

5.00 

.1890  -2 

76.920 

10.00 

.2356  -4 

102.32 

5.10 

.1683  -2 

77.841 

10.10 

.2205  -4 

102.59 

5.20 

.1501  -2 

78.733 

10.20 

.2065  -4 

102.85 

5.30 

.1341  -2 

79.597 

10.30 

.1934  -4 

103.11 

5.40 

.1200  -2 

80.434 

10.40 

.1813  -4 

103.36 

5.50 

.1075  -2 

81.245 

10.50 

.1701  -4 

103.61 

5.60 

.9643  -3 

82.032 

10,60 

.1596  -4 

103.86 

5.70 

.8663  -3 

S2.795 

10.70 

.1499  -4 

104.10 

5.80 

.7794  -3 

83.537 

10.80 

.1408  -4 

104.33 

5.90 

.7021  -3 

84.257 

10.90 

.1324  -4 

104.57 

9352  x 

10"1  -  .9352 

-1 

TABLE  18.  (Continued) 


M 

P/Pt 

V 

M 

P/Pt 

V 

11.00 

.1245  -4 

104.80 

16.30 

.8565  -6 

113.00 

11.10 

.1171  -4 

105.02 

16.40 

.8213  -6 

113.11 

11.20 

.1103  -4 

105.24 

16.50 

.7876  -6 

113.21 

11.30 

.1039  -4 

105.46 

16.60 

.7556  -6 

113.31 

11.40 

.9788  -5 

105.67 

16.70 

.7250  -6 

113.41 

11.50 

.9228  -5 

105.88 

16.80 

.6959  -6 

113.51 

11.60 

.8704  -5 

106.09 

16.90 

.6680  -6 

113.61 

11.70 

.8215  -5 

106.29 

17.00 

.6415  -6 

113.71 

11.80 

.7755  -5 

106.49 

17.10 

.6161  -6 

113.81 

11.90 

.73  25  -5 

106.69 

17.20 

.5918  -6 

113.90 

12.00 

.6922  -5 

106.88 

17.30 

.5687  -6 

114.00 

12.10 

.6544  -5 

107.07 

17.40 

.5465  -6 

114.09 

12.20 

.6189  -5 

107.26 

17.50 

.5254  -6 

114.18 

12.30 

.5857  -5 

107.44 

17.60 

.5052  -6 

114.27 

12.40 

.5544  -5 

107.62 

17.70 

.4859  -6 

114.36 

12.50 

.5250  -5 

107.80 

17.80 

.4674  -6 

114.45 

12.60 

.4973  -5 

107.98 

17.90 

.4496  -6 

114.54 

12.70 

.4714  -5 

108.15 

18.00 

.4328  -6 

114.63 

12.80 

.4469  -5 

108.32 

18.10 

.4165  -6 

114.72 

12.90 

.4239  -5 

108.49 

18.20 

.4010  -6 

114.80 

13.00 

.4023  -5 

108.65 

18.30 

.3861  -6 

114.89 

13.10 

.3818  -5 

108.82 

18.40 

.3718  -6 

114.97 

13.20 

.3626  -5 

108.97 

18.50 

.3582  -6 

115.05 

13.30 

.3444  -5 

109 . 13 

18.60 

.3452  -6 

115.13 

13.40 

.3273  -5 

109.29 

18.70 

.3326  >6 

115.21 

13.50 

.3111  -5 

109.44 

18.80 

.3206  -6 

115.29 

13.60 

.2958  -5 

109 . 59 

18.90 

.3090  -6 

115.38 

13.70 

.2814  -5 

109.75 

19.00 

.2980  -6 

115.45 

13.80 

.2678  -5 

109.89 

19.10 

.2874  -6 

115.53 

13.90 

.2550  -5 

110.04 

19.20 

.2772  -6 

115.61 

14.00 

.2428  -5 

110.18 

19.30 

.2674  -6 

115.68 

14.10 

.2313  -5 

110.32 

19.40 

.2581  -6 

115.76 

14.20 

.2204  -5 

110.46 

19.50 

.2491  -6 

115.83 

14.30 

.2100  -5 

110.60 

19.60 

.2404  -6 

115.91 

14.40 

.2003  -5 

110.74 

19.70 

.2321  -6 

115.98 

14.50 

.1910  -5 

110.87 

19.80 

.2241  -6 

116.05 

14.60 

.1823  -5 

111.00 

19.90 

.2165  -6 

116.13 

14.70 

.1739  -5 

111.13 

20.00 

.2091  -6 

116.20 

14.80 

.1660  -5 

111.26 

14.90 

.1586  -5 

111.38 

15.00 

.1515  -5 

111.51 

15.10 

.1447  -5 

111.63 

15.20 

.1383  -5 

111.76 

15.30 

.1323  -5 

111.88 

15.40 

.1265  -5 

112.00 

15.50 

.1210  -5 

112.11 

15.60 

.1158  -5 

112.23 

15.70 

.1108  -5 

112.34 

15.80 

.1061  -5 

112.45 

15.90 

.1016  -5 

112.57 

16.00 

.9731  -6 

112.68 

16.10 

.9323  -6 

112.79 

16.20 

.8936  -6 

112.89 
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For  a  -  0°,  sin  6  ■  sin  0. 

The  value  for  C  on  the  desired  meridian,  Or,  may  be  determined 
^max 

either  experimentally  or  theoretically.  The  NASA*1  ,sa  zero  and  small  angle 

M 

of  attack  cone  tables  can  be  used  for  the  starting  values  on  a  cone  fitted  to 
the  nose  vertex  if  experimental  data  are  lacking;  the  procedure  is  as  follows: 
The  general  expression  for  pressure  coefficient  is: 


c  -  (  JL  -  1 )  and  f  - 

'  %\Po  I  Po  p,  Po 


where  p,  refers  to  conditions  at  zero  angle  of  attack  and  the  quantity  p, /pQ 
may  be  obtained  directly  from  the  zero  angle  of  attack  cone  tables**  ,  The  ratio 
of  p/pj  which  is  the  ratio  of  static  pressure  at  angle  of  attack  to  static 
pressure  at  zero  angle  of  attack  can  be  calculated  using  the  theory  of  Stone*3 
wherein  the  velocity,  pressure,  and  density  are  expanded  in  the  following  series 
and  higher  order  terms  in  a  are  neglected: 

H*  -  M,*  -  av£  cos  * 

P  “  Pi  -  ops  cos  * 

p  ■  p,  -  a  pa  cos  * 

where  the  flow  quantities  lif ,  Pj  ,  and  p:  refer  to  conditions  at  jsero  angle  of 

attack  and  14  ,  ps  ,  and  pa  are  the  flow  quantities  related  to  the  effect  of  angle 

of  attack.  Reference  23  provides  solutions  for  the  above  equations  yielding, 


9 - 

The  necessary  tables  required  in  the  calculational  procedures  of  this  sub¬ 
section  have  been  reproduced  from  the  references  and  Included  in  Appendix  1. 


among  others,  the  following  expression: 


P/Pi 


(r+1)  -  (r-1)  M** 
(r+D  -  (r-D  mT* 


a 


e 


(S./R)  cos  iji 


References  21  and  22  have  provided  tabulated  values  of  Mi*,  t£ , 

and  Sa/R-for  use  in  determining  p/p5  and  subsequently,  Cp.  Since  conditions  at 

th^  nose  vertex  were  used  to  calculate  the  pressure  coefficient,  this  value  of 

C  becomes  C  from  whence  C  ,  Equation  (7),  may  be  determined.  Starting 
p  PN  pmax 

values  of  C  must  be  calculated  for  each  value  of  *  at  the  desired  angle  of 
pmax 

attack.  Once  this  is  done,  the  Generalized  Newtonian  theory,  Equation  (6),  is 
used  to  calculate  the  values  of  Cp  along  each  meridian  up  to  the  point  where  the 
shock  expansion  method  is  used  to  extend  the  impact  theory;  for  the  cylindrical 
afterbody,  Equation  (5a)  or  (5b)  is  used  for  the  pressure  coefficient  distribut¬ 
ion. 


The  local  Mach  number  at  the  nose  vertex  (where  M  -  V  may  be 
determined  from  the  relation*4: 


II*  -  2  M*a 

(r+D  -  (y-l)  M** 

On  the  leeward  side  of  the  body,  certain  portions  of  the  surface 

will  lie  in  the  "aerodynamic  shadow"  and  the  Newtonian  theory  cannot  predict  the 

pressures  in  this  region.  This  shadow  area  on  the  body,  shown  pictorlally  here, 


can  be  determined  quite  easily  along  any  leeward  meridian  by  setting  sin  6  -  0 
and  solving  for  0.  For  any  meridian  in  general,  the  0  at  which  the  shadowed 
area  starts  is: 

0  -  tan"1  (tan  a  cos  *) 
for  *  -  0°  0  -  a 
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Thus  the  shadow  area  begins  along  the  ♦  ■  0°  meridian  at  the  point  where  the 
body  surface  angle  is  equal  to  tiie  angle  of  attack. 

On  all  meridians,  the  matching  point  (where  the  shock-expansion 
method  is  used  to  extend  the  impact  theory)  is  obtained  directly  from  plots 
presented  in  prior  sections  of  this  report,  e.g.,  Pig.  6,  page  16.  When  the 
matching  point  lies  within  the  shadowed  area,  the  shock-expansion  method  should 
be  started  at  the  x/t  value  where  sin  6-0. 

A  numerical  example  is  presented  in  Appendix  1  wherein  the  pressure 
coefficient  distribution  has  been  calculated  for  a  tangent  ogive-cylindrical 
afterbody  combination  at  MQ  -  2.0  and  1/ d  -  3  at  an  angle  of  attack  of  5  degrees. 
Tables  are  included  which  provide  the  necessary  parameters  for  determining  the 
pressure  and  local  Mach  number  distribution.  For  the  numerical  example,  start¬ 
ing  values  of  C  were  determined  for  seven  meridians  beginning  at  ♦  -  0*  and 
PN 

proceeding  in  30  degree  increments  up  to  *  -  180°.  The  pressure  distribution 
along  the  +  -  180°  meridian  is  the  only  one  for  which  the  Cp  calculation  is  shown 
since  the  procedure  is  similar  for  each  meridian.  The  calculr.tion  of  the  cylin¬ 
drical  afterbody  pressure  distribution  using  Equation  (5b)  has  also  been  in¬ 
cluded  for  the  *  -  180°  meridian. 

The  pressure  distributions  calculated  by  the  present  method  for  the 
numerical  example  of  Appendix  1  are  presented  in  Figure  23  where  experimental 
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Fig.  23  PRESSURE  COEFFICIENT  DISTRIBUTION  FOR  A  TANGENT 
OGIVE-CYLINDRICAL  AFTERBODY  COMBINATION  AT  a  -  5°, 

*  -  0°,  30°,  60°,  90J,  120°,  150°,  and  180° 


values'8  have  been  included  for  comparison  purposes  with  the  subsequent  indi¬ 
cation  of  good  agreement. 

2.  Bodies  of  Revolution  at  Large  Angles  of  Attack 

The  calculation  of  pressure  distributions  for  angles  of  attack 

■X-  ^ 

greater  than  about  S  degrees  requires  somewhat  more  time  since  the  starting 

values  of  C  must  be  determined  from  the  kit**'*7’*8  cone  tables  (or  comparable 
PN 

information)  unless  experimental  data  are  available.  Once  the  starting  values 
have  been  obtained,  the  procedure  is  similar  to  that  used  for  the  small  angles 
of  attack. 

Since  the  MIT  cone  tables  are  being  used  to  obtain  starting  values 

at  the  nose  vertex,  the  symbols  and  nomenclature  adopted  by  this  widely  used 

reference  will  not  be  altered  to  conform  to  the  nomenclature  of  this  report. 

Any  attempt  to  redefine  or  re-reference  the  parameters  would  inevitably  tend  to 

compound  the  existing  complexity  of  the  equations.  The  symbols  and  nomenclature 

from  the  MIT  cone  tables  which  are  used  in  the  present  study  are  defined  on  page 

4  of  this  report.  In  addition,  the  parameters  required  for  calculating  C  have 

PN 

been  picked  from  the  MIT  cone  tables  and  are  tabulated  in  Appendix  2  for  ready 
use.  If  any  intermediate  values  of  the  parameters  are  required,  they  can  be 
obtained  much  more  readily  when  nomenclature  consistent  with  the  reference  source 
is  used. 

The  theory  of  Stone*3 is  used  once  again  to  determine  the  flow 
parameters  which  have  been  expanded  in  the  following  series;  for  large  angles 
of  attack,  the  higher  order  terms  in  a  cannot  be  neglected,  i.e., 

P/p"  ■  1  +  a  A,  cos  *  +  a*  (A*  +  As  cos  2t) 

c/o  “  1  +  a  B,  cos  t  +  a*  (B,  +  B3  cos  2*) 

where  p/p  and  p/p"  are  the  pressure  and  density  on  the  cone  surface  at  angle  of 
attack  divided  by  the  corresponding  values  at  zero  angle  of  attack,  A,  and  B, 
specify  the  first  order  effects  of  a,  and  A|  ,  As,  B, ,  and  Bs  represent  the  second 
order  effects  of  a.  Proceeding  further: 

A,  -  -ri/p 

^  -  Po^  +  -5"  fr  +  2^  COt  ®s 

As"  Ps/P  +  - - —  cot  ©„ 

*  a*  2p  8 
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Fig.  24  PRESSURE  COEFFICIENT  DISTRIBUTION  FOR  A  TANGENT 
OGIVE-CYLINDRICAL  AFTERBODY  COMBINATION  AT 
a  -  10°,  *  -  0°,  30°,  60®,  90s,  120®,  150",  and  180® 
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